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ABSTRACT 


The  kinetic  mechanisms  responsible  for  the  excitation  of  the  first 
positive  and  first  negative  emission  of  nitrogen  have  been  investigated  in 
a  re-examination  of  previously  reported  shock-tube  measurements  of  the  non 
equilibrium  radiation  for  these  systems.  The  rate  coefficients  of  the 

colllslonal  quenching  reactions: 

.00 


and 


N.(A  3I+)  +  N(AS)  — N„  (X  1E+)  +  N(4S) 
/  u  2  g 


t(N2) 


N*(B  2E+)  +  K,(X  1E+) - — >N?(X  2E+)  or  N^(A  2H  )  +  N,(X  1£+) 

zu^g  2  g  2u2g 

were  found  to  be  given  by  the  empirical  expressions: 

.00  ,  .  ,_-3  -2.2'j  3  -1 

k_2  ■  5.1  x  10  T  cm  sec 

and 

i  (No)  .  0  in-2  -2.33  3  -1 

k  *  ■  1,9  x  10  T  cm  sec  , 

q  ’ 

respectively,  over  the  approximate  temperature  range  6000°  -  14000°K. 
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1.  INTRODUCTION 

The  characteristic  nonequilibrium  radiation  profiles  which  ere 

observed  in  the  relaxation  region  behind  shock  waves  in  nitrogen  in  the 

o  o 

nominal  wavelength  intervals  6000-12000  A  and  3500-5000  A  are  generally 

3  3  + 

attributed  to  the  first  positive  (B  II  A  E  )  and  first  negative 

8 

2  +  2  +  + 

(B  E^  -*•  X  1^)  band  systems  of  N2  and  N2,  respectively.  The  characteris¬ 
tic  features  of  this  nonequilibrium  emission  have  been  investigated  by 
several  authors.^  ^  Qualitatively,  the  radiation  rises  rapidly  to  a 
peak  immediately  after  the  shock  wave,  the  peak  intensity  being  far  in 
excess  of  that  which  would  correspond  to  complete  thermochemical  equili¬ 
brium,  and  then  decays  in  an  exponential-like  fashion  to  the  ultimate 

equilibrium  level.  The  emitting  states  of  the  two  band  systems,  l.e., 

3  +  2  + 

N2(B  T!  )  and  ^(B  E^) ,  are  separated  by  more  than  10  e.V,  and  thus  it 

is  surprising  that  the  times  required  to  reach  the  peak  intensity  are 

quite  comparable  over  a  wide  range  of  temperature,^  ^  typically  6000°-15000°K. 

This  observation  suggests  common  rate-limiting  steps  in  the  kinetic 

mechanisms  which  are  responsible  for  populating  both  of  the  emitting  states 

N- (B  3n  )  and  N*(B  2E+). 

^8  /  u 

In  an  attempt  to  examine  the  role  of  atomic  nitrogen  in  shock-wave 
excitation  mechanisms,  Wray  ^ ^  fired  shock  waves  into  a  nitrogen  test 
gas  which  was  already  partially  dissociated  by  means  of  a  pulsed  electrode¬ 
less  bischarge.  Unfortunately,  the  initial  atom  concentrations  in  the 
test  gjs  at  the  time  when  the  shock  waves  passed  the  observation  station 
were  not  accurately  determined.  The  primary  source  of  the  Inaccuracy 
appears  to  be  that  the  low  temperature  recombination  rate  coefficient  which 


-3- 


Wray  used  Co  estimate  the  atom  concentration  at  the  time  of  shock  arrival 
was  roughly  a  factor  of  ten  greater  than  the  more  recent  and  generally 
accepted  value.  However,  on  the  basis  of  Wray's  measurements, 

it  is  possible  to  conclude  that  nitrogen  atoms  are  very  much  more  effec¬ 
tive  in  promoting  collisional  excitation  of  nitrogen  molecules  to  both  the 
N2(B  ng)  and  N2(B  states  than  are  the  ground  state  molecules  h'2(  Eg). 

Recent  improvements  in  our  understanding  of  the  gas  phase  dissocia¬ 
tion  and  recombination  kinetics  of  nitrogen,  both  experimentally^8 
and  theoretically, prompted  us  to  re-examine  the  role  of  atomic  nitro¬ 
gen  in  the  excitation  mechanisms  which  give  rise  to  the  first  positive 
and  first  negative  emission  profiles  observed  behind  shock  waves.  In 
this  investigation  we  have  used  the  measurements  of  the  times  to  reach 
peak  intensity  and  the  radiation  profile  shapes  which  have  been  published 

in  the  literature, together  with  some  unpublished  measurements,  due 
(12) 

to  Wray,  of  the  times  to  reach  peak  intensity  and  the  absolute  emission 

intensities  recorded  in  experiments  where  the  test  gas  was  not  predissociated. 

2.  EXPERIMENTAL  CONSIDERATIONS 

2.1  vibrational  Relaxation 
(12) 

Figure  1  shows  a  plot  of  Wray's  measurements  of  the  times,  T^, 

to  reach  peak  intensify  for  both  the  N2(l+)  and  N2(l-)  band  systems  as 

a  function  of  temperature.  All  of  the  experiments  were  conducted  using 

undiluted  nitrogen  at  an  initial  pressure  of  1  torr.  The  full  line  in 

Figure  1  shows  the  temperature  variation  of  the  characteristic  vibrational 

relaxation  time  for  ground  state  N2  as  measured  by  Millikan  and  White^18^ 

(19) 

and  Appleton.  It  is  clear  from  this  comparison  that  the  nitrogen  gas 
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was  vibrationally  relaxed  well  before  the  N2(I+)  and  N^l-)  radiation 
intensities  had  reached  significant  levels;  thus  the  initial  post-shock 
temperatures  were  calculated  using  the  Rankine-Hugoniot  equations,  and  the 
assumption  that  the  translational,  rotational,  and  vibrational  energy 
modes  of  the  N2  were  fully  equilibrated.  However,  with  increasing  dis¬ 
tance  behind  the  shock  front,  the  temperature  and  pressure  will  vary  due 
to  both  the  dissociation  of  N2  molecules  and  to  the  development  of  the 
cold  wall  boundary  layer  which  acts  as  a  mass  sink  for  the  shocked  test 
gas. 

2.2  Boundary  Layer  Effects 

The  effects  of  the  shock-tube  boundary  layer  on  the  properties  of 

f20-26) 

the  shocked  test  gas  have  been  investigated  by  numerous  authors.  ' 

For  this  reason  we  shall  not  dwell  on  the  subject  here.  Mirel 's^25'26^ 
treatment,  which  allows  corrections  to  the  ideal  shock-tube  flow  proper¬ 
ties  due  to  boundary  layer  development,  is  the  most  widely  quoted.  The 
primary  effect  of  the  boundary  layer  growth  in  kinetic  studies  performed 
in  small  bore  shock  tubes  at  low  initial  pressures  is  to  decrease  the  time 
of  flight  of  a  fluid  element  which  passes  through  the  shock  wave  at  a 
fixed  position  upstream  of  the  observation  station  below  that  calculated 
on  the  basis  of  ideal  shock  tube  theory.  The  corrections  to  the  ideal 
flow  properties  are  usually  evaluated  by  assuming  that  at  the  observation 
station  the  shock  and  contact  surface  are  travelling  at  the  same  speed; 
i.e,,  the  mass  leakage  in  the  boundary  layer  which  passes  the  contact 
surface  is  equal  to  the  mass  flux  through  the  shock  front.  This  condi¬ 
tion  was  approximately  satisfied  in  Wray's  experiments.  The  mass  continuity 
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equation  for  the  inviscid  flow  external  to  the  boundary  layer  may  then 
be  written  In  the  form 

pu  *  P2u2(l  -  (x/Am)1^2)  (2.1) 

where  P2^2  *  p^Ug,  is  the  mass  flux  entering  the  shock  wave,  x  is  the 
distance  measured  downstream  from  the  shock  wave  in  shock-fixed  coordi¬ 
nates,  and  is  the  maximum  separation  distance  between  the  shock  and 
the  contact  surface.  Since  it  is  assured  that  the  shock  waves  travel  at 
the  uniform  velocity  Ug,  then 

x  -  Dt.  (2.2) 

8  L 

where  tT  Is  the  laboratory  observation  time* 

In  order  to  analyze  the  chemically  reacting  flows  obtained  behind 
normal  shock  waves  generated  in  a  shock  tube,  we  have  developed  a  numeri¬ 
cal  computer  program  which  uses  a  fourth-order  Runga  Kutta  Integration 
technique  to  solve  the  differential  forms  of  the  flow  conservation  equa¬ 
tions  together  with  the  appropriate  chemical  rate  equations.  When 
viewed  in  shock-fixed  coordinates,  the  flow  is  steady  and  is  treated  as 
being  quasi-one-dimensional  to  account  for  the  boundary  layer  mass  loss 
as  described  by  Equation  2.1.  The  independent  variable  used  in  the  analysis 
is  x  which  is  related  to  the  actual  particle  flight  time  t^  by  the  equation 
dx/dtf  -  u  .  (2.3) 

For  a  more  detailed  discussion  of  the  calculation  method,  see  the  Appendix. 

2.3  Emission  Intensity  Calculation 
Having  determined  the  temperature  density  and  species  concentra¬ 


tion  histories  in  the  relaxation  region,  the  spectral  band  intensities 
of  the  radiation  were  calculated  using  the  sine  ire d  rotational  line  model 


— 6— 


(4) 

described  by  Keck,  Camm,  Kivel,  and  Wentink.  Thus  the  spectral  Intensity 
1^  is  given  as 

IA  -  2hc2(nrQ)  f  [N]  <<j)>  X“5(hc/kT)  g ' 'Qr ' 7g'Qr ’C^’  (2.4) 

2  2  _  2 
where  r  -  e  /me  is  the  classical  electron  radiur>,  f  -  |R(r)/ea  |  /3R  X 
o  o 

_  2 

is  the  absorption  oscillator  strength  of  the  band  system,  and  |R(r)/eaQ| 
is  the  electronic  transition  moment  dependent  on  the  internuclear  separa¬ 
tion  r.  The  quantity  <J>  is  a  dimensionless  number  of  order  unity  which 
takes  into  account  the  details  of  the  vibration-rotational  spectrum,  and 
the  averaged  quantity  is  defined  as 

O)-  X6  /2  (X')“6  RdX'//2  RdX‘  (2.5) 

A1  A1 

where  R  is  the  experimentally  determined  resolution  function  of  the  opti¬ 
cal  instruments  used  in  the  experiments.^11^  [N'J  is  the  upper  state 
concentration,  and  g',  Q^',  and  Q  '  are  the  corresponding  electronic 
degeneracy,  rotational,  and  vibrational  partition  functions,  respectively, 
which  are  evaluated  by  assuming  that  the  translational,  rotational,  and 
vibrational  temperatures  are  equal  (the  double  primes  identify  the  corres¬ 
ponding  absorbing  state  quantities). 

In  the  calculations  which  were  carried  out  for  the  purpose  of  compari- 

n  2) 

son  with  Wray's  absolute  intensity  measurements,  values  for  the  elec- 

(27) 

tronic  transition  moments  were  those  given  by  Wurster  for  the  N^(l+) 

i  _  2  (28) 

system,  i.e.,  |R(r)/eao|  ■  0,096,  and  by  Buttrey  and  McChesney  for 

,  2 

the  N2(l-)  system,  i.e.,  |R(r)/eao|  -  0.45, 

3.  KINETIC  MODELS  AND  COMPUTATIONS 

3,1  Reaction  Mechanism  -  ( 1+) 

The  kinetic  scheme  which  we  have  found  to  best  describe  the  shock- 
tube  measurements  of  the  N2(l+)  system  is  the  following: 
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1  +  k(M) 

N2(X  hj  +  M  2n(4s)  +  M 

k(N) 

n9(x  1z+)  +  n(4s)  n9(a  3z*)  +  n(4s) 


,00 

-2 


u 


(1) 

(2) 


N_  (A  3E+)  +  H 
2  u 


k(M) 

*■+•»  4 

2N(  S)  +  M 


(3) 


k(M) 

N_(A  3E+)  +  M  v--  ->  N,(B  3n  )  +  m 

u  k<M)  8 

-4 


(4) 


X“1 

n„(b  3n  ) — ^->n„(a  3r+)  +  hv 

2  g  2  u 


(5) 


Apart  from  our  inclusion  of  reaction  (3),  the  above  set  of  reactions  is 


the  same  as  that  previously  proposed  by  Wray. 


(ID 


The  overall  rate  of 


the  dissociation  of  ground  state  N2  behind  shock  waves  has  been  measured 
over  the  temperature  range  8000°-15000°K  by  Appleton,  Steinberg,  and 


Liquornik 


(16) 


for  those  cases  where  the  collision  partner  M  is  either 

Ar,  N2,  or  N,  Shui,  Appleton,  and  Keck,  ^)  using  the  modified  phase- 

space  theory  of  reaction  rates,  have  calculated  the  separate  dissociation 

(M)  (m) 

rate  coefficients  k+^  and  k+^  for  the  case  where  the  collision  partner 

3  + 

M  is  an  argon  atom.  By  assuming  that  the  N2(A  £  )  state  was  in  local 
thermodynamic  equilibrium  with  the  ground  state,  they  were  able  to  match 
both  the  absolute  magnitude  and  temperature  dependence  of  the  measured 
dissociation  rate  coefficients  quite  well  using  what  appears  to  be  a 
fairly  realistic  form  for  the  two-body  interatomic  potential  On 

the  basis  of  the  results  to  be  described  later,  it  appears  that  for  the 
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dissoclatlon  of  N2  diluted  in  an  argon  heat  bath  and  at  the  temperatures 
of  the  shock-tube  experiments,  the  assumption  of  local  thermodynamic 
equilibrium  is  invalid.  We  shall  also  show  that  for  the  dissociation  of 
pure  N2  at  temperatures  greater  than  about  S000  K,  the  A  E^  state  con¬ 
centration  is  well  below  that  calculated  on  the  basis  of  local  thermody¬ 
namic  equilibrium  during  the  rise  time  of  N2(l+)  radiation,  and  thus  the 
dissociation  of  both  undiluted  N2  and  of  diluted  N2  at  high  temperatures 
behind  shock  waves  proceeds  primarily  via  the  ground  state,  i.e.,  reaction 
(1).  We  have  therefore  used  the  experimentally  derived  dissociation  rate 

coefficient  k^  for  k^;  but  have  additionally  assumed  that  the  ratio 

(M)  (M) 

k+3  /k+1  (M  -  i!2,  Ar,  h)  is  the  same  as  that  given  by  the  phase-space 

theory  calculations  for  M  »  Ar.^^  ihe  ultimate  justification  for  this 

assumption  is  provided  by  the  results  contained  below;  however,  even  on 

theoretical  grounds  we  should  anticipate  this  procedure  to  be  approximately 
* 

correct. 

In  accord  with  Wray,^11^  we  have  assumed  reaction  (4)  to  be  suffi¬ 
ciently  fast  by  comparison  with  any  of  the  other  reactions  which  serve 


3  +  3 

to  populate  or  depopulate  either  of  the  N2(A  Ey)  or  N2(B  11^)  states 
that  these  two  states  may  be  assumed  to  be  in  local  thermodynamic  equi¬ 
librium  with  one  another  throughout  the  entire  relaxation  region.  This 

assumption  was  subsequently  justified  by  a  numerical  calculation  in  which 

(N) 

we  used  our  derived  estimate  for  k+2  ,  the  measured  radiative  lifetime  of 
3  (39)  -6 

the  B  H  state  (T,  ■  7  x  10  sec),  and  an  estimate  of  the  collxslonal 
g  5 


The  phase-space  theory  assumes  that  the  three-body  interaction  potential 
is  given  by  the  sum  of  two  two-body  potentials  Vj;2  and  where  V^2  Is  the 

ground  state  (X  1E+) aolecular  potential  for  the  calculation  of  k^O ,  and  the 
first  excited  state  (A  3£+)  molecular  potential  for  the  calculation  of  k+g'. 
Vjftj  is  simply  the  Interatomic  potential  between  an  N(^S)  atom  and  the  third 
body  M;  it  is  independent  of  the  ultimate  state  of  the  combined  atoms. 
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(M)  “XI  3  —1 

quenching  rate  coefficient  k_^  ■  6.2  x  10  cm  sec  obtained  at 

T  *  300°K.  The  rate  of  the  dissociation  reaction 

N0(B  3JI  )  +  M  N(4S)  +  N(2D)  +  M 

2  g 

was  shown  theoretically^2^  to  be  negligible  by  comparison  with  reac¬ 
tion  (3)  due  to  its  Increased  endothermlcity. 

3.2  Calculations  and  Comparison  with  Experiment.  ^(1+) 

Using  the  above  values  for  k+^  ,  ,  and  (see  Table  I  for  a 

summary  of  the  individual  rate  coefficients  used)  and  the  assumption  that 

3  +  3 

the  A  £u  and  B  11^  states  were  maintained  in  local  thermodynamic  equi¬ 
librium,  we  were  able  to  deduce  the  N-atom  excitation  rate  coefficient, 

(N)  (N) 

k+2  "  Kg2  k_2  i  by  matching  our  computed  ^(1+)  emission  profiles  with 

those  observed  experimentally.  The  matching  procedure  used  was  as  follows: 

The  rise  of  the  emission  intensity  behind  the  shock  wave  was  calculated  to 

the  peak  using  a  guessed  estimate  for  the  value  of  the  rate  coefficient 

k^^  •  The  computed  maximum  intensity  was  then  compared  with  Wray's^2^ 

absolute  measured  value.  By  iteration  on  the  value  of  the  rate  coeffi- 
(N) 

clent  k_2  ,  the  calculated  peak  Intensity  was  brought  into  agreement  with 
the  experimental  value.  This  matching  procedure  was  employed  at  several 
shock  speeds  which  spanned  the  full  experimental  range  of  conditions.  The 
empirical  rate  coefficients  thus  obtained  were  correlated  by  the  following 
expression: 


,  (N)  c  .  in-3  -2,23  3  -1 

k_2  ■  5,1  x  10  T  cm  sec 


(3.1) 


for  the  temperature  range  6000°-14000°K.  Figure  2  shows  the  comparison 

.  (12) 

of  our  calculated  peak  Intensities  with  Wray  s  measurements. 
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(N) 

In  Figure  3  we  have  compared  the  values  of  k_2  given  by  Equation  3.1 

with  those  obtained  by  Wray^^  for  the  same  temperature  range,  and  with 

(31) 

the  room  temperature  rat**»  determined  by  Young  and  St.  John  and  by 

(32)  (N) 

Meyer,  Setser  and  Stedman,  It  is  apparent  that  our  estimate  of  k  2 

Is  more  than  an  order  of  magnitude  greater  than  Wray's  estimate.  We 
shall  discuss  the  comparison  shown  in  Figure  3  later  in  this  report. 

Additional  tests  of  the  kinetic  model,  reactions  (1)  -  (5),  and  of 
the  rate  coefficients  are  provided  by  comparisons  of  the  theoretically 
calculated  times-to-peak  intensity  and  the  intensity  profile  shapes  with 
the  corresponding  measurements;  this  information  was  not  used  in  the 
matching  procedure  described  above.  Figure  4  presents  a  compilation  of 
tlme-to-peak  Intensity  data  for  the  first  positive  emission  obtained  In 
shock-tube  experiments  using  undiluted  N2/Ar,^  and  N2/Ne^^ 

mixtures.  It  is  apparent  that  the  calculated  values  of  T  ^  for  pure 
nitrogen  agree  very  well  with  the  experimental  data.  Similar  good  agree¬ 
ment  Is  observed  for  the  N2/Ar  mixtures,  although  the  recorded  time  was 

not  the  actual  time  to  peak  but  rather  a  characteristic  time,  t'  ,  as 

pic 

shown  In  Figure  5.  Since  we  anticipate  that  the  dissociation  rate  coeffi¬ 
cients,  k^®^  and  k^e\  for  neon  as  the  collision  partner  do  not  differ 

(Ar)  (Ar) 

greatly  from  the  argon  rate  coefficients,  k+^  and  k^  *  the  actual 
times  to  peak  calculated  for  the  N2/Ar  mixtures  have  also  been  plotted 
for  comparison  with  the  measurements  made  in  N2/Ne  mixtures.  Again  the 
agreement  between  theory  and  experiment  appears  to  be  satisfactory. 

A  representative  emission  intensity  profile  shape  is  presented  In 
Figure  5  for  an  observation  made  behind  a  shock  wave  in  pure  N2>  The 
shaded  region  represents  the  relative  noise  level  of  the  oscilloscope 
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trace,  and  the  full  line  represents  the  theoretically  calculated  profile. 
The  similarity  between  the  theoretical  and  experimental  profiles  which 
extends  well  Into  the  region  where  thermochemical  equilibrium  is  approached, 
further  substantiates  the  kinetic  model  and  the  rate  coefficients  used. 

The  calculated  concentration  histories  shown  in  Figure  6a  and  b  for 
a  shock  wave  in  undiluted  Nj  and  in  a  10  per  cent  N^/Ar  mixture,  respec¬ 
tively,  are  helpful  to  our  understanding  of  the  radiation  overshoots  of 
the  first  positive  emission.  The  pure  nitrogen  case,  Figure  6a,  corres¬ 
ponds  to  the  radiation  intensity  profile  of  Figure  5.  Due  to  the  shock- 
tube  boundary  layer  development,  the  ground-state  molecule  concentration 
is  observed  to  increase  slightly  with  time.  The  atom  concentration 
remains  small,  reaching  only  10  per  cent  of  the  ground-state  molecule 
concentration  at  times  well  after  the  peak  intensity  is  achieved,  although, 
of  course,  for  stronger  shock  waves  the  degree  of  dissociation  increases 
rapidly,  becoming  about  25  per  cent  at  the  peak  intensity  for  a  post¬ 
shock  temperature  of  about  14000°K, 

3  + 

The  A  Ey  state  concentration  at  the  intensity  peak  is  about  one- 

quarter  of  the  local  thermodynamic  equilibrium  value  which  is  represented 

by  the  dashed  curve  in  Figure  6a.  This  difference  increases  with  increas- 

3  + 

ing  temperature  so  that  the  ^(A  Eu)  concentration  at  peak  intensity  is 

less  than  one-tenth  the  equilibrium  value  at  T£  ■  14000°K,  whereas,  at 

m  6000°K,  it  is  greater  than  one-third  of  the  equilibrium  ccncentra- 

3  + 

tion.  The  reason  for  this  relative  behavior  of  the  A  E  concentration 

u 

can  be  understood  in  terns  of  the  expression  which  describes  its  steady- 

3  + 

state  variation.  The  steady  state  ^(A  E^)  concentration  is  achieved 

at  about  the  time  of  the  intensity  maximum,  and  is  given  by: 
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(N*)  -  (N,)  (tO/lkj^)  (f,2)  +  k|N)  (N)  +  k(N)  (N)]  (3.2) 

*  3  +x 

where  N2  ■  N2(A  Zy) .  At  high  temperatures  where  the  atom  concentration 
is  large  at  the  peak  and  the  rate  kj_^  is  significantly  greater  than  k^\ 
the  second  term  of  the  denominator  dominates  since  the  atom  is  much  more 
efficient  than  the  nitrogen  molecule  as  a  collision  partner  in  the  dissocia¬ 
tion  reactions.  Thus  Equation  3.2  may  be  approximated  by 

(N*)  -  k£j}  (N2)/k|^  .  (3.3) 

At  low  temperatures  (T2  <  6000°K)  where,  due  to  the  exponential  behavior 
of  k^j\  the  term  k^^  (N)  at  the  peak  intensity  is  much  greater  than  k|^^ 
(N,? J ,  the  steady-state  approximation  reduces  to  the  local  thermodynamic 
equilibrium  condition: 

* 

(N*)  -  Ke2  (N2)  .  (3.4) 

The  experimental  shock  tube  observations  of  the  N2(l+)  emission  only 
approach  these  two  extremes  of  the  temperature  range  which  satisfy  either 
Equations  3.3  or  3,4.  Over  much  of  the  experimental  temperature  range, 
the  term  k^2^  (N2)  is  largest,  but  because  the  additional  terms  of  the 
denominator  are  not  wholly  nv'.gliglble,  the  steady-state  concentration  is 
only  approximated  by 

(N*)  -  kj2}  (N)/k^2>  .  (3.5) 

A 

By  tacitly  assuming  that  during  the  intensity  rise  the  N  and  ^ 
concentrations  are  populated  according  to  the  simplified  rate  equations 

d(N)/dt  -  2k|j2)  (N2)2 


(3.6) 
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d(N*)/dt  «  k^  (N2)  (N) 


(3.7) 


and  that  the  relaxation  process  Is  isothermal  and  the  amount  of  dissocia¬ 
tion  is  so  small  that  the  net  rates  of  reactions  (1)  and  (2)  remain 
constant,  we  obtain  by  integration 


(N)  -  2k^2>  (n2)2  t 


(3.8) 


and 


(N*)  -  kj^  (N2,3  t2  . 


(3.9) 

We  must  remember  that  for  pure  nitrogen,  particularly  at  high  temperatures, 

the  assumptions  of  constant  H2  concentration  and  temperature  are  valid 

only  for  times  considerably  shorter  than  the  peak  time.  However,  by 

* 

approximating  the  population  of  the  N2  state  by  Equation  3.9  until  the 
steady-state  concentration  is  reached,  the  important  features  of 
become  apparent.  At  the  extreme  high  temperatures  where  the  steady-state 
concentration  is  given  by  Equation  3.3,  we  find 


<V  <N2)rl  '  <k«2>  k«>1/2  • 


(3.10) 


At  the  low  temperature  extreme  where  the  assumptions  of  this  simplified 
analysis  are  more  closely  satisfied,  the  local  equailibrium  concentra¬ 
tion  is  reached  at  the  peak,  giving 


(TpkOyr1  *  (k|^)  k^jV/2 . 


(3.11) 


The  steady-state  concentration  over  the  intermediate  experimental  tempera¬ 
ture  range  is  given  approximately  by  Equation  3.5;  thus 

(t  (N  ))”^  *\>  k^2^ 

'  pk'  2"  +3 


(3.12) 
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Since  the  assumptions  leading  to  Equations  3.8  and  3.9  are  only 

approximations  at  best,  we  may  only  reasonably  expect  the  observed  time 

to  peak  intensity  to  be  correlated  by  a  plot  of  log  [  (t  ^  (N2)  )”*] 

versus  T  \  as  suggested  by  the  form  of  Equations  3.10-3.12.  Indeed, 

this  expectation  appears  to  be  fully  realized  by  the  results  shown  in 

Figure  4.  However,  as  we  mentioned  above,  over  most  of  the  experimental 

temperature  range,  i.e.,  T  >  6000°K,  the  N2(A  2E*)  concentration  is  much 

smaller  than  the  corresponding  equilibrium  concentration,  and  thus  the 

calculated  times  to  peak  are  primarily  determined  by  the  dissociation 

reactions  (1)  and  (3)  and  are  much  less  dependent  on  reaction  (2). 

Therefore,  the  good  agreement  between  the  calculated  and  measured  times 

to  peak  as  shown  in  Figure  4  must  reflect,  the  essential  accuracy  of  the 

(M) 

values  which  we  have  used  for  the  dissociation  rate  coefficients  k^ 
and  for  M  •  N2  and  N. 

3  + 

Figure  6b  shows  that  the  N2(A  state  concentration  for  a  dilute 

N2/Ar  mixture  remains  well  below  the  corresponding  equilibrium  concentra¬ 
tion  even  beyond  the  intensity  maximum.  A  similar  analysis  as  that  out- 

3  + 

lined  above  for  the  steady-state  concentration  of  the  A  state  shows 
that  in  the  limit  of  a  dilute  mixture,  it  is  approximately  given  by 

[N*3  -  kj2}  [N2]  [N]/k^r)  [ArJ  (3.13) 

over  the  experimental  shock-tube  temperature  range.  Using  equations 
similar  to  Equations  3,8  and  3.9,  we  obtain 

(T^lAr])"1  *  kJAr)  .  (3.14) 

The  same  observation  which  was  made  for  the  case  of  undiluted  N2  is  again 


valid;  i.e.,  t  e  good  agreement  between  the  calculated  and  measured 
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tlmes-to-peak  lntesnlty  shown  In  Figure  4  for  dilute  mixtures  s  -bstantl- 

(Ar) 

ate  the  essential  accuracy  of  the  value  for  k+^  which  we  have  used  In 
our  calculations. 

(N' 

The  accuracy  of  the  rate  coefficient  k _2'  which  we  have  determined 

here  is,  of  course,  dependent  on  the  validity  of  the  assumed  kinetic 

scheme,  the  accuracy  of  the  other  rate  coefficients  which  we  have  used, 

and  the  reliability  of  the  flow  model  which  has  been  used  to  take  account 

of  the  shock-tube  boundary  layer  effects.  We  have  carried  out  exploratory 

calculations  in  which  we  altered  the  values  of  the  rates  of  reactions  (1) 

and  (3)  by  factors  of  two;  we  also  performed  the  calculations  without 

making  the  boundary  layer  calculations.  In  this  way  we  sought  to  establish 

(N) 

the  sensitivity  of  our  derived  values  for  k_2  to  the  various  assumptions. 

(N) 

On  this  basis  we  conclude  that  the  expression  for  k_2  (Equation  3.1)  is 
reliable  to  within  a  factor  of  about  two  over  the  temperature  range  6000°- 
14000°K.  Although  there  is  no  experimental  evidence  on  which  to  base 
estimates  of  k^?p  within  the  temperature  range  300°-600Q°K,  it  is  clear 
from  the  results  shown  in  Figure  3  that  a  straight-line  extrapolation  would 
correlate  both  the  room  temperature  measurements  and  our  shock-tube  estimates 
quite  well. 

Before  leaving  the  subject  of  the  N2(l+)  emission,  we  shall  briefly 

comment  on  two  other  reactions  which  have  been  suggested  as  being  primary 

3  + 

sources  of  the  excited  state  N„(A  E  ).  The  first 

2  u 

2N,(X  1I+)  -*•  N_ (A  V)  +  N.(X  1I+)  (2a) 

2  g  2  u  2  8 

is  spin  forbidden,  and  although  the  N2  ground  state  is  by  far  tne  most 


abundant  species  immediately  behind  the  shock  waves,  we  were  unable  to 
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determine  a  rate  coefficient  for  this  reaction  by  our  matching  procedures 

which  enabled  us  to  correlate  both  the  peak-intensity  measurement  and  the 

(33) 

time-to-peak  intensity.  In  addition,  it  is  to  be  noted  that  Noxon 

concluded  from  observations  made  in  high  pressure  nitrogen  afterglow  experi- 

9 

ments  that  more  than  10  collisions  with  ground  state  molecules  are  required 

3  + 

to  effect  electronic  de-excitation  of  the  A  Z  state.  Thus  we  believe 

u 

that  the  above  excitation  reaction  is  too  slow  to  be  of  importance  in  the 
shock  wave  experiments, 

Smekhov  and  Losev v  '  have  suggested  an  alternative  mechanism  for  the 
excitation  of  N^(A  3E*) 


N_ (X  1Z+)  +  e" 
2  g 


===^N_(A  3I+)  +  e" 

e  u 


where  the  electrons  are  generated  by  the  associative  ionization  reaction 

(e) 

to  be  discussed  in  the  following  section.  Although  the  rate  k  2  nay  be 
(N) 

much  larger  than  k_2  ,  it  is  unlikely  that  the  difference  will  be  as  great 
as  the  difference  between  the  N  and  e  concentrations  except  at  very  high 
temperatures.  Using  a  recent  estimate  of  the  rate  coefficient  of  reaction 
(2b),(3A)  l.e. , 


.(e)  +  o\  in“5  -0.8  -80300/T  3  -1 

k^2  »  (5  -  2)  x  10  T  e  cm  sec 


and  our  cal  :uiated  values  for  the  concentrations  and  temperatures  at  the 
N2(l+)  emission  peaks,  we  conclude  that  only  as  T2  approaches  15000°K 
does  the  contribution  of  reaction  (2b)  reach  one-tenth  of  the  rate  of 
production  of  N2(A  3L^)  by  N2  +  N  collisions. 


3.3  Reaction  Mechanism  -  N+d-) 

Th'i  kinetic  scheme  which  we  have  found  to  best  describe  fie  shock- 

tube  measurements  of  the  N^d-)  system  follows  directly  from  that  proposed 

(3) 

by  Hammerling,  Tea re,  and  Klvel: 


N(4S)  +  N(4S)  ^====*N*(X  2E*)  +  e" 

k-6 


2  g 


N_ (X  XE+)  +  N+(3P)  ^====iN*(X  2E+)  +  N(4S) 
2  g  u  2  8 

-7 


(M) 


+  2  + 

»T(x  i  ) 

2  g 


4<a  2nu) 


■v - 

■ 

k(M) 

-8 

k«0 

+9 

— 

k(M) 

-9 

k(M) 

+10^ 

V 


2V  V 


Nl(B  "E  )  +  M 


n*(x2e+)+m^ - , 

8  k(M)  2 

-10 


(6) 

(7) 

(8) 

(9) 

(10) 


+  2  +  '  1 1 1  +  2  + 

N^B  Eg)  +  hV  (11) 

(35) 

Dunn  and  Lord!  have  recently  estimated  the  associative  ioniza¬ 
tion  rate  coefficient  k  ^  over  the  approximate  temperature  range  3500°- 
7200°K  from  measurements  of  the  electron  density  decay  rate  obtained  In 
a  shock-tube  wind  tunnel  nozzle.  We  have  assumed  that  the  expression 
which  they  give  for  k_^  can  be  extended  to  the  higher  temperatures  which 


are  of  Interest  to  us  here. 
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Char  ge-exchange  cross  sections  of  the  type  associated  with  reaction 
(7)  are  known  to  be  quite  large;  however,  only  crude  estimates  of  the 
rate  coefficients  k+^  -  k_^  have  been  made;  for  example,  Dunn  and 

Lordi^^  suggest  k  ^  *  1,3  x  10^  ^  cm^  sec  An  estimate  of 

the  rate  coefficient  for  reaction  (8),  i.e.,  k^g2^  •  10  ^  cm^  sec  \ 

f  £_  \ 

has  been  given  by  Bennett  and  Dalby  at  T  ■  300CK;  we  anticiprte 

that  k  g  would  not  exhibit  a  significantly  stronger  temperature  depend- 

(N) 

ence  than  was  observed  for  k  ^  »  see  Figure  Thus  we  carried  out 

trial  calculations  based  on  the  estimates  of  k  ^  and  k  g  given  above 

and  concluded  that  reactions  (7)  and  (8)  would  achieve  local  equilibrium 

well  within  the  time  required  for  the  N+d-)  radiation  to  reach  its 

intensity  maximum  over  the  experimental  shock-tube  temperature  ran/,e. 

+  2 

To  proceed  further,  we  assumed  that  not  only  were  the  ^(A  IO 
+  3 

and  N  (  P)  states  in  local  equilibrium  with  the  ground  states  of  the 

molecular  ion,  but  in  addition,  we  also  assumed  that  the  concentration 
2 

of  the  ^(B  E^)  state  was  given  by  the  steady-state  approximation  at 

all  times;  i.e., 


[N+(B  V)J 


u  (M) 

+9 


IN* (A  2nu)] 


[M]  +  kff  [N+(X  V)]  [M] 


k^}  [M]  +  kfj*  (M}  +  t'J 


•  Kel0‘B2«  !  ^ 


k(M)  [M] 


+  Tn 


(3.16) 


2  +  Q  /qn 

where  the  radiative  life-time  of  the  B  E  state  is  T,..  ■  6.58  x  10  sec, 

u  11 

K  m[Nt(x  2^+)J  is  the  equilibrium  concentration  of  the  B  2E+  state  relative 
eiu  L  g  u 

to  the  ground  state  of  the  ion,  and  k^^  ■  (k^^  +  k^g)  is  the  net  colli- 

2  + 

slonal  quenching  rate  of  the  B  E^  state  which  we  sought  to  determine. 
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3.4  Calculations  and  Comparison  with  Experiment:  NjCl-) 
Radiation  intensity  profiles  for  the  I^Cl-)  system  were  first  com¬ 
puted  for  a  range  of  conditions  which  spanned  Wray's  experimental  shock- 
tube  measurements.  The  dissociation  kinetics  of  the  were  taken  to 
be  the  same  as  those  described  previously,  ionization  was  assumed  to 

proceed  via  reaction  (6),  reactions  (7)  and  (8)  were  assumed  to  be  in 

+  2  + 

local  equilibrium,  and  the  N?(B  E^)  state  concentration  was  assumed  to 

be  given  by  the  steady-state  condition,  Equation  3.16.  By  iteration  on 
(M) 

the  value  of  at  each  temperature,  it  was  a  straightforward  matter 

to  match  the  calculated  peak  intensity  with  the  absolute  peak  intensity 

(12) 

measurements  given  by  Wray.  In  this  way  we  deduced  the  following 

empirical  expression  for  the  net  collisional  quenching  rate  coefficient 
2  + 

for  the  B  E  state: 


k(N2> 

q 


1.9  x  i(f2  t"2,33 


cm  sec 


-1 


(3.17) 


The  calculated  peak  intensities  are  compared  with  Wray's  measurements  in 

Figure  7,  and  a  plot  of  k^N2^  (Equation  3.17)  is  shown  in  Figure  3  together 

with  the  estimates  of  the  low  temperature  quenching  rate  coefficients 

determined  by  Brocklehurst ^32^  and  Davidson  and  O'Neil^3^  at  T  *  300°K. 

Figure  8  shows  the  comparison  between  measurements  of  the  time  required 

to  reach  the  intensity  maximum  and  our  calculated  estimates  of  these  times. 

+  2  + 

Since  we  have  assumed  that  the  ^(B  Ey)  state  concentration  was  given  by 
the  steady-state  approximation,  the  calculated  values  of  T  ^  are  virtually 
independent  of  the  rates  of  reactions  (8)  and  (9).  Rather,  the  good  agree¬ 
ment  between  the  experimental  and  theoretical  values  of  T  ^  shown  in 
Figure  8  reflect  the  essential  accuracy  of  the  dissociation  rates  and  the 
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high  temperature  extrapolation  of  the  rate  for  the  associative  ionization 
reaction,  reaction  (6),  used  in  the  calculations.  Also  shown  in  Figure  8 
is  the  comparison  between  measurements  of  the  characteristic  time  T by 
Smekhov  and  Losev ^  (open  circles)  for  10  per  cent  N^/Ar  mixtures  and 
computed  values  (dashed  line)  for  which  we  assumed  k^r^  ■  Again 

the  comparison  is  seen  to  be  quite  good. 

It  is  to  be  noted  that  the  comparisons  shown  in  Figure  8  are  effected 
by  plotting  log  {t  [M])^^}  *  versus  T  *  rather  than  log  * 

versus  T  *  as  in  Figure  4  for  the  ^(1+)  system.  We  shall  now  examine 
the  reason  for  this  method  of  comparison. 

The  species  concentration  profiles  illustrated  in  Figures  6a  and  b 
show  that  the  dominant  ionic  species  up  to  the  time  of  the  intensity  maxi¬ 
mum  is  the  ground  state  molecular  ion.  Furthermore,  our  calculations 
indicated  that  the  peak  Intensity  was  achieved  at  about  the  time  at  which 
local  equilibrium  was  approximately  established  for  reaction  (6);  i.e., 

(3.18) 


[N+(X  V))  =  K^2  [N]  . 


By  assuming  that  during  the  major  portion  of  the  period  the  ^  con¬ 

centration  did  not  change  significantly  and  that  the  conditions  behind 
the  shock  wave  are  Isothermal,  we  obtain 


-  kj”)[N2]  [M] 

and 

d[N^]  2 

ir  •  k«  i*” 

for  t  <  T  .  ,  and  thus 

pk 


(3.19) 


(3.20) 
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[N*]  -  4k^)2  k+6[N2]2  [Mj  t3/3  .  (3.21) 

By  further  assuming  that  Equation  3.21  is  valid  up  to  the  point  at  which 
local  equilibrium  is  approximately  established  for  reaction  (6) ,  we  then 
obtain 

{tpk([N2]  [M])1/2}-1  *  (kj^0  k_6  K^2)1/2  .  (3.22) 

This  approximate  result  explains  why  T  ^  is  primarily  determined  by  the 
dissociation  kinetics  of  N2  and  the  rate  of  the  associative  ionization 
reaction,  as  previously  suggested,  and  the  reason  for  our  choice  of  the 
form  for  the  ordinate  in  Figure  8. 

Figure  9  shows  a  comparison  between  a  representative  emission 
intensity  profile  (shaded  region)  for  an  observation  behind  a  shock  wave 
in  an  N2/Ar  mixture  with  the  theoretically  calculated  profile  given  by 
the  full  line.  The  similarity  between  the  theoretical  and  experimental 
profiles  further  substantiates  the  kinetic  model  and  rate  coefficients 
used. 


As  we  observed  for  the  case  of  the  N2(l+)  system,  the  reliability 

+  2  + 

of  our  estimate  for  the  collisional  quenching  rate  of  the  N2(B  1^)  state 

is  dependent  on  the  validity  of  the  model  reaction  scheme  and  on  the 

accuracy  of  the  other  rate  coefficients  used  in  the  calculations.  For 

this  model  reaction  scheme,  we  have  adjusted  various  of  the  rates  over 

their  probable  range  of  uncertainty  to  assess  their  relative  importance 

in  the  overall  calculation,  and  on  this  basis  we  suggest  that  our  estimate 
(M) 

of  kq  is  probably  accurate  to  within  a  factor  of  2  or  3.  Again,  we 

point  out  that  a  linear  extrapolation  from  our  high  temperature  estimates 
(M) 

of  k^  does  appear  to  correlate  with  the  low  temperature  measurements 
shown  in  Figure  3. 
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4.  CONCLUSIONS 


A  computer  program  has  been  developed  which  allows  calculations  to 
be  made  of  chemically  reacting  flows  obtained  behind  normal  shock  waves 
produced  in  non-ideal  shock  tubes.  With  the  aid  of  this  program,  non¬ 
equilibrium  radiation  intensity  profiles  for  the  N2(l+)  and  N+Q-)  band 

systems  have  been  calculated  and  compared  with  previously  published 

(12) 

and  unpublished  measurements.  The  model  reaction  schemes  used  In 

the  calculations  followed  directly  from  those  proposed  in  earlier  studies, 

and  were  found  to  yield  good  agreement  with  measurements  of  the  absolute 

maximum  intensity,  time-to-peak  Intensity,  and  intensity  profile  shapes. 

3  + 

Rate  coefficients  for  the  colllslonal  quenching  of  the  N^(A  £^)  and 
+  2  + 

N2(B  E^)  states  via  the  reactions: 

k(N) 

n2(a  V)  +  n(4s)  “2-  >n2(x  V)  +  n(as) 

and 

k<N2> 

N,(B  h+)  +  N„ - 9 - >N,(A  2n-  )  or  N,(X  V)  +  N 

2  U  2  2  u  2  g  2 

(12) 

were  deduced  by  matching  the  computed  peak  Intensities  with  Wray's 
measurements  and  were  found  to  be  correlated  by  the  empirical  expressions 


,  (N)  c  ,  -3  -2.23  3  -1 

k_2  -  5.1  x  10  T  cm  sec 

and 

.  (No)  i  n  ,»-2  -2.33  3  -1 

k  2  ■  1.9  x  10  T  cm  sec 

q 

over  the  temperature  range  6000°-140n°K.  The  two  rates  so  determined 


are  observed  to  be  of  comparable  magnitude, and  both  have  a  similar  strong 
negative  temperature  dependence.  The  reasons  for  this  similarity  are  not 
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clear,  although  It  Is  to  be  noted  from  Figure  3  that  these  high  tempera¬ 
ture  estimates  of  the  rate  coefficients  extrapolate  quite  well  to  the 
corresponding  low  temperature  measurements. 

The  good  agreement  between  the  computed  and  measured  time-to-peak 
intensity  and  the  general  profile  shapes  for  both  the  ^(1+)  and  N^Cl-) 
systems  was  observed  to  result  from  the  choice  of  the  rate  coefficients 

which  describe  the  dissociation  kinetics  of  nitrogen^^*^"^  and  of  the 

(35) 

ionization  mechanism.  In  particular,  it  appears  that  Shul,  Appleton, 

and  Keck's estimate  of  the  ratio  k^/k^,  is  substantiated,  although 

3  + 

it  is  apparent  that  their  assumption  that  the  ^(A  Z^)  and  ground  states 

are  in  local  equilibrium  is  invalid  for  dissociation  behind  shock  waves 

at  temperatures  greater  than  about  6000°K.  Since  the  contribution  of 
3  + 

the  A  Z(_  state  to  the  net  dissociation  rate  is  smaller  than  they  suggested, 

it  explains  the  lack  of  an  observed  induction  time  in  the  nitrogen  dissocia- 
(39) 

tion  measurements.  Finally,  as  a  consequence  of  the  comparisons  pre¬ 

sented  here  for  the  N+d-)  system,  it  appears  that  Dunn  and  Lordi's^^ 
expression  for  the  rate  coefficient,  k+^,  of  the  associative-ionization 
reaction  can  be  extended  up  to  14000°K. 

5.  APPENDIX 

Nonequilibrium  Shock-Tube  Program 

The  flow  properties  behind  normal  shock  waves  in  a  reacting  gas  are 
determined  by  integration  of  the  differential  forms  of  the  flow  conserva¬ 
tion  equations  for  the  steady  quasl-one-dlmensional  flow  as  observed  in 
a  shock-fixed  coordinate  system: 

djO  P2U2  1  _  £  du 

dx  2u  ,  .  vl/2  u  dx 

(xE  ) 
m 


(5.1) 
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dp  „  du 
■7*1  ■  -  pu  — — 
dx  dx 


dh  du  n 

dx  dx 


(5.2) 

(5.3) 


The  effects  of  shock-tube  boundary  layer  development  on  the  inviscid  flow 
are  taken  into  account  by  the  "mass  sink"  term  which  appears  on  the  right- 
hand  side  of  Equation  5.1.  The  subscript  (2)  refers  to  the  conditions 
Immediately  behind  the  shock  front. 

By  defining  Y^  ■  N^/p,  where  is  the  number  concentration  of  species 

1,  h ,  as  the  enthalpy  per  molecule  of  species  1  so  that  h  **  E  y  h  ,  and 
1  1  1  1 

Cp^  ■  dh^/dT,  Equation  5.3  becomes 


v  u  ^1  .  dT  v  _  du  _ 

L  h.  -7--  +  -r-  L  y  Cp .  +  u  7—  *  0  . 

^  1  dx  dx  ^  1  *1  dx 


(5.4) 


For  a  perfect  gas 


RP 


where  the  molecular  weight  M  nay  be  expressed  as 
N 


M 


£  Y, 


(5.5) 


(5.6) 


By  differentiating  Equation  5.5  with  respect  to  x,  we  find 

dT  m  _  P_  M_  y  M_  dP  _  PM_  d£ 

dx  Rp  N  .  dx  Rp  dx  D  2  dx 

O  1  Kp 


where 


dx  pu  dt 

and  dN^/dt  is  g^ven  by  the  usual  kinetic  rate  equation. 

5.9  may  now  be  combined  !:c  yield 


(5.8) 

Equations  5.1  - 
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•“idT  d».  T(^ic^ 

pu  +  »oP«  (l  \  Cpl)  dt  >  -  2(x,  ,1/2  (1  ;  (xl  ,1/2, 

m  in 


T  E  Yi  CPi  2 

u  +  - (1  -  £2-> 


(5.9) 


Equations  2.3,  5.2,  5.8,  and  5.9  fotm  the  basis  for  a  numerical  Integra¬ 
tion  which  yields  the  complete  nonequilibrium  flow  properties  in  the 
reaction  zone  behind  the  shock  front. 


Reaction  Rate  Coefficient*  Reference 


i 


TABLE  I.  Reaction  Rate  Coefficients  Used  in  Calculation  of  Radiation  Profiles 
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FIGURE  CAPTIONS 

(12) 

Fig.  1  Comparison  of  measured  times-to-peak  intensity  due  to  Wray 
(P^  ■  1  Torr)  and  the  characteristic  vibrational  relaxation 
time,  tv.<18-19)  •  ,  N2(l+);  O,  N+U-);  full  line,  Ty. 

Fig.  2  Comparison  of  measured  and  calculated  peak  intensities  of  the 
N2(l+)  emission.  •,  Wray;  full  lines,  this  work. 

Fig.  3  Comparison  of  collisional  quenching  rates,  k^:  #,  Young 
and  St.  John;  A,  Meyer,  Setser,  and  Stedman;^2^  full 
lines,  this  work  and  Wray^x^  as  indicated.  k^N2^;  O, 

q 

Brocklehurst;  A,  Davidson  and  O'Neil;  broken  line, 
this  work. 

Fig.  4  Comparison  of  measured  and  calculated  time-to-peak  intensity  of 

(12) 

the  first  positive  system.  Wray:  #,  100%  N2,  P1  «  1  Torri 

Wurster  and  Marrone:^  A,  100%  N2,  Pj  -  1-^5  Torr;  Keck,  Camm, 
Kivel,  and  Wentink:  ,  100%  N2>  P^  -  1  Torr;^^,  P^  **  3  Torr; 

®  ,  P^  ■  10  Torr;  Smekhov  and  Losev:  ^  O,  10%  N2  +  90%  Ar, 

P^  »  2-10  Torr;  Marrone,  Wurster,  and  Stratton:  #,  10% 

N2  +  50%  Ne,  Px  -  2  Torr;  72  N2  +  932  Ne»  Pi  “  2  Torr;  A,  32 

N2  +  97%  Ne,  Px  -  2  Torr;  O  ,  3%  N£  +  97%  Ne,  -  4  Torr;  V,  0.5% 

N2  +  99.5%  Ne,  Px  -  4  Torr;  () ,  3%  N  +  97%  Ne,  -  6  Torr;  ▼ , 

0.5%  N2  +  99.5%  Ne,  Px  -  6  Torr. 

_ ,  calculated  T^,  100%  N2,  P^  ■  1  Torr;  -  -  - - -  , 

calculated  T^,  10%  N2  +  90%  Ar,  P^  ■  6  Torr; - , 

calculated  T  .  ,  10%  N_  +  90%  Ar,  P,  ■  6  Torr. 
pk  2  *  1 
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Fig.  5  Comparison  of  theoretically  calculated  N^(l+)  emission  profile 
(full  line)  with  experimental  profile  of  Wray.  100% 

■  4.56  mm/ysec,  **  1  Torr,  T2  ■  8210  °K. 

Fig.  6  Computed  concentration  histories.  _  ,  calculated; 

_  ,  local  thermodynamic  equilibrium;  _  , 

temperature. 

a)  Wray^*^  profile:  100%  Ug  ■  4.56  mm/ysec,  ■  1  Torr, 

T2  -  8210  °K. 

b)  Smekhov  and  Losev^  profile:  10%  N_  +  90%  Ar,  U  ■  3.35  mm/ysec, 

A  8 

Px  -  6  Torr,  T2  -  9550  °K. 

Fig.  7  Comparison  of  measured  and  calculated  peak  intensities  of  the 

N+d-)  emission.  •,  Wray; ^ 12 ^  full  line,  this  work. 

Fig.  8  Comparison  of  measured  and  calculated  time-to-peak  intensity  of 

the  first  negative  system.  Wray:  #,  100%  N2,  P^  ■  1  Torr; 

Smekhov  and  Losev: ^  O,  10%  N2  +  90%  Ar,  P^  -  1-10  Torr. 

Full  line:  x  .  ,  100%  N»;  broken  line:  x  ,  ,  10%  N„  +  90%  Ar. 
pk  2  pk  2 

Fig.  9  Comparison  of  theoretically  calculated  N^d-)  emission  profile 
(full  line)  with  experimental  profile  (shaded  region)  due  to 
Smekhov  and  Losev: ^  10%  N2  +  90%  Ar,  Ug  ■  3.35  mm/ysec, 

Pj  -  6  Torr,  T2  -  9550  °K. 


«*  -I 


TIME  TO  PEAK  INTENSITY  r„k  (sec) 


6 


8 


FIGURE  5 


CONCENTRATION  (molecules/cm3) 


TEMPERATURE  (°KxlO* 


CONCENTRATION  (molecules /cm 


TEMPERATURE  (°K  x  10 


PEAK  INTENSITY  ( watts  /cm3  •  ster  •  jj. ) 


FIGURE  9 


